Gap-anisotropic model for the narrow-gap Kondo insulators 
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A theory is presented which accounts for the dynamical generation of a hybridization gap with 
nodes in the Kondo insulating materials CeNiSn and CeRhSb. We show that Hunds interactions 
acting on virtual 4/ 2 configurations of the cerium ion can act to dynamically select the shape of 
the cerium ion by generating a Weiss field which couples to the shape of the ion. In low symmetry 
crystals where the external crystal fields are negligible, this process selects a nodal Kondo semimetal 
state as the lowest energy configuration. 

Pacs numbers: 



Kondo insulators share in common with the Mott insu- 
lators a gap which is driven by interaction effects. 1 ' 2 Un- 
like Mott insulators, they undergo a smooth cross-over 
into the insulating state, where a tiny charge and spin 
gap develops. These materials are generally regarded as 
a special class of heavy fermion system, where a lattice 
Kondo effect between the localized spins and conduction 
electrons forms a highly renormalized band-insulator. 3 ' 4 
The smallest gap Kondo insulators, CeNiSn and 
CeRhSb, do not naturally fit into this scheme: they 
appear to develop gapless excitations. Early measure- 
ments showed a drastic increase of the electrical resistiv- 
ity below 6K, 5 but very pure samples of CeNiSn display 
metallic behavior. 6 NMR measurements are consistent 
with an electronic state with a "v-shaped" component to 
the density of states. 7 These results, together with other 
transport properties 8-11 point to the formation of a new 
kind of semi-metal with an anisotropic hybridization gap. 

Ikeda and Miyake 12 (IM) recently proposed that the 
Kondo insulating ground-state of these materials devel- 
ops in a crystal field state with an axially symmetric hy- 
bridization potential that vanishes along a single crystal 
axis. This picture accounts for the v-shaped density of 
states, and provides an appealing way to understand the 
anisotropic transport at low temperatures, but it leaves a 
number of puzzling questions. In CeNiSn and CeRhSb, 
the Cerium ions are located at sites of minimal mono- 
clinic symmetry, where the low-lying f-state is a Kramers 
doublet 

|±) = hi ± 1/2) + 6 2 | ± 5/2) + 6 3 | T 3/2) (1) 

where b = (61,62,63) could point anywhere on the unit 
sphere, depending on details of the monoclinic crystal 
field. The IM model corresponds to three symmetry- 
related points in the space of crystal field ground states, 

6 = <^ ^+T~> 4 ' I> (2) 
1(0,0,1) [> 

where a node develops along the x, y or z axis respec- 
tively. What mechanism selects this special semi-metal 



out of the manifold of gapped Kondo insulators? Neu- 
tron scattering results show no crystal field satellites in 
the dynamical spin susceptibility of CeNiSn, 13 suggesting 
the the crystal electric fields are quenched: is the selec- 
tion of the nodal semi-metal then a many body effect? 14 

In this letter, we propose that this selection mechanism 
is driven by Hund's interactions amongst f-electrons in 
the Cerium ions. Hund's interactions play an important 
role in multi f-electron ions. 15 In the Kondo semi-metal, 
the Cerium ions are in a nominal 4/ 1 state, but undergo 
valence fluctuations into f° and f 2 configurations. We 
show that the memory effect of the Hund's interactions 
in the f 2 state induces a kind of Weiss field which couples 
to the shape of the Cerium ion. When this field adjusts 
to minimize the Hund's interaction energy, the nodal IM 
state is selected. 

To develop our model, we classify each single-particle 
f-configuration by a "shape" (a=l,2,3) and a pseudo-spin 
quantum number(ev = ±1), where 

/+ ± |0> = |±l/2>, 
fl ± \0) = \±5/2), 

/] ± |0)^|t3/2) (3) 
There are eight multipole operators 

[T]° = flAlJca, (« = 1,8) (4) 

which describe the shape of the Cerium ion, where the 
A a matrices are the eight traceless SU(3) generators, nor- 
malized so that Tr[A a A b ] = S ab , We shall describe the 
low energy physics by an Anderson model H — H + Hj, 
where 

H = H c + ]T V[cUj)faAj) + H.c], (5) 

jaa 

and H c = ^ kcr e^c^ c ka describe a spin-1/2 conduction 
band hybridized with a lattice of localized f-states. The 
operator 

cL(j) = (N.yi Xy^y^kjcL, (6) 

k,cr 
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creates a conduction electron hi a I — 3, j — b/2 Wannicr 
state at site j with shape-spin quantum numbers (a, a), 
N s is the number of sites and 

y: a (k) = r 3 " lJ ^(k)(^ °, 3 mj - a\\mj) (7) 

defines the form-factors, in terms of spherical harmon- 
ics and the Clebsh-Gordon coefficients of the j = 5/2 
/ 1 -state 16 , where mj = mj(a,a) maps the spin-shape 
quantum numbers to to the original azimuthal quantum 
number of the f-scattering channel. Following previous 
authors, 17 we regard if as a low energy Hamiltonian, so 
that hybridization strength V is a renormalized quantity, 
that takes into account the high energy valence and spin 
fluctuations. 
The term 

H f =^[E f n f (j) + ^(n f (j) - l) 2 - §r 2 ] (8) 

j 

describes the residual low-energy interactions amongst 
the f-electrons: the second term is a Coulomb interac- 
tion term. The third term is a Hunds interaction which 
favors 4/ 2 states with maximal total angular momentum. 
In an isotropic environment, this interaction would take 
the form — | J 2 , where J is the total angular momentum 
operator, but in a crystalline environment, it takes on a 
reduced symmetry which we model in simplified form by 
— |r 2 . In general the Hund's interaction is only invari- 
ant under discrete rotations so that fluctuations into the 
f 2 state enable the system to sample the crystal symme- 
try even when the conventional crystal field splittings are 
absent. 

Suppose the crystal electric field term were un- 
quenched, so that H — > H — ^2 a ■ Tj. The shape of 
the Cerium ion (Ij ) = T is determined by the condition 
that the energy is stationary with respect to variations 
in T, 

N- 1 6(H )/5T = a + gT. (9) 

The second term is a feedback or "Weiss" contribution 
to the crystalline electric field, created by fluctuations 
into the 4/ 2 state. Generally, the induced field T will 
follow the crystalline electric fields a, but in situations 
where the valence and spin fluctuations are rapid enough 
to quench the external crystal electric field, 13 then a = 0, 
and the Weiss field becomes free to explore phase space 
to minimize the total energy. In such a situation, the 
shape of the Cerium ion is determined by the interac- 
tions, rather than the local conditions around each ion. 

To explore this process, we carry out a Hubbard- 
Stratonovich decoupling of the interactions, 

Hf(j) - f] [(E f + A,)l + Aj • A] f s + E [Xj,Aj] (10) 
where 

Eo [\ j ,A j } = &-§--\ j ), (11) 



Here A a (j) <~ — gT a (j) is a dynamical Weiss field, fj 
denotes the spinor fj = f a a{j)- Note that in the path 
integral, the fluctuating part of Xj, associated with the 
suppression of charge fluctuations, is imaginary. We now 
seek a mean-field solution where the Weiss field Xj = A 
and Aj = A, and E(Xj,Aj) = E Q . Such an expectation 
value does not break the crystal symmetry. However, the 
selected crystal field matrix A-A must adjust to minimize 
the total energy. Supposing we diagonalize this matrix, 
writing A • A = UA W, where Aq = diag(Ai, A 2 , A 3 ), 

and Ai > A 2 > A 3 . In the basis, f au (j) = Ul b f b(7 (j), 
the crystal field is diagonal. In practice, the strength of 
the Hund's interaction g is so large that the excitation 
energies Ai i2 — A 3 substantially exceed the Kondo tem- 
perature. In this case, the mean-field Hamiltonian must 
be projected into the subspace of the lowest eigenvalue. 
In the hybridization, we therefore replace 

c)fi = c]Ufj -+ b a [cUj)U(j)}, (12) 

where f a {j) = fs<y (dropping the superfluous index "3") 
describes the lowest Kramers doublet and b a = C/ a3 . To 
satisfy the constraint (n/) = 1, the energy of the lowest 
Kramers doublet must be zero, i.e. Ef + A + A 3 = 0. We 
then arrive at the mean-field Hamiltonian 

H* = H c + V^[^ ra (k)4/ak + H.c] + N S E (13) 

k 

where 4> aa (k) = ^ a b a y% a (k) is the dynamically gener- 
ated form-factor of the hybridization. 16 The transformed 
hybridization is no longer rotationally invariant: all infor- 
mation about the anisotropic wavefunction of the Cerium 
ion is now encoded in the vector b. 

The quasiparticle energies associated with this Hamil- 
tonian are 

£ k ± = £k /2± v /( ek /2) 2 + V k 2 (14) 

Here, the hybridization can be written in the con- 
venient form V£ = V /2 $j ) (k) where ^(k) = 

(V 2 ) E Q . CT I E Q b a yi a (k)| 2 contains all the details of the 
gap anisotropy. The ground-state energy is then the sum 
of the energies of the filled lower band 

E 9 = - 2 E \/( £ k/2) 2 + V* + N S E (15) 

k 

Now both A and A 3 are fixed independently of the direc- 
tion of b, so that E does not depend on b. To see this, 
write the eigenvalues of the traceless crystal field matrix 
as Ai i2 = "^A ± S, A 3 = — -^A. Since the upper two 
crystal field states are empty, stationarity w.r.t to S re- 
quires <5 = 0. Since A 3 couples directly to the f-charge, 

we obtain dE g /dA = -y^n/) + (A/g) = 0, so that 
A = J\g. Thus both A = -A 3 - E f and A Q arc fixed 
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Global Minima 

FIG. 1. Contour plot of the ground-state energy in 
mean-field theory as a function of the two first components 
of the unit vector b (the third one is taken as positive). The 
darkest regions correspond to lowest values of the free-energy. 
Arrows point to the three global and three local minima. 



independently of b. The selection of the crystal field con- 
figuration is thus entirely determined by minimizing the 
kinetic energy of the electrons. To examine the depen- 
dence of the mean-field on b, we replace the momentum 
sum in (15) by an energy and angular integral, 



(16) 



where N(0) and 2D are respectively, the density of states 
and band- width of the conduction band. Completing the 
integral, noting that the angular average ($g(k)) = 1, 
we find that the shift in the ground-state energy per site 
due to the hybridization is 



AE g = 2N(0)V 2 



ln[ 



eD 2i 



F[b] 



where 



F[b] 
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$ g (k)ln[$ g (k)] 



(17) 



(18) 



The weak logarithmic divergence inside F(h) favors 
states with nodes. Fig. 1. shows a contour plot of the 
mean-field free-energy as a function of the two first com- 
ponents of b. There are three global minima and three 
local minima with a slightly higher frcc-cncrgy. The state 
where b = z, plus two symmetry equivalents, corresponds 
to the IM state and has the lowest free-energy. The 
IM state is axially symmetric, with a hybridization node 
along the z, y or x axis. But the theory also identifies 
a new locally stable state where b = (0, -\/5/4, \/TT/4), 
plus its two symmetry equivalents. This state is almost 
octahedral. Like the IM state, the hybridization drops 



exactly to zero along the z axis. But, in marked dif- 
ference with the IM state, it almost vanishes along the 
(1, 1, 0) and (1, — 1, 0) directions in the basal plane. 

The relative stability of the IM and the octahedral 
state will, in general be dependent on details of our 
model, such as the detailed conduction electron band- 
structure. For this reason, both possibilities should 
be considered as candidates for the nodal semi-metallic 
states of CeNiSn and CeRhSb. The inset in Fig (2) 
shows the density of states predicted by these two possi- 
bilities. Although both are gapless, the v-shaped pseudo- 
gap of the quasi-octahedral state is far more pronounced 
than in the axial state, and is closer in character to 
the observed tunneling density of states. 18 A more di- 
rect probes of the anisotropy is provided by the ther- 
mal conductivity 19 which, unlike the resistivity, does not 
show a strong sample dependence in these compounds. 
To compute and compare the theoretical thermal conduc- 
tivity with experiments, we compute the thermal current 
correlator 20 



1 

2T 



— - / dbjuj' 



df\N(cj) 
duj J T(lu) 



(19) 



where / is the Fermi function, T(lu) is the quasiparticle 
scattering rate and 



(20) 



describes the quasiparticle velocity distribution where 
= VgEg and Eg is given by equation (14). For our 
calculation, we have considered quasiparticle scattering 
off a small, but finite density of unitarilly scattering im- 
purities or "Kondo holes". 21 We use a self-consistent T- 
matrix approximation, following the lines of earlier cal- 
culations except for one key difference. In these calcu- 
lations, which depend critically on the anisotropy, it is 
essential to include the momentum dependence of the 
hybridization potential in the evaluation of the quasipar- 
ticle current. Previous calculations 12 underestimated the 
anisotropy by neglecting these contributions. 20 

The single node in the IM state leads to a pronounced 
enhancement of the low-temperature thermal conductiv- 
ity along the nodal z axis. By contrast, in the quasi- 
octahedral state the distribution of minima in the gap 
give rise to a modest enhancement of the thermal conduc- 
tivity in the basal plane. Experimental measurements 19 
tend to favor the latter scenario, showing an enhancement 
in thermal conductivity that is much more pronounced 
in k x than in k z or K y . 

Three aspects of our theory deserve more extensive ex- 
amination. Nodal gap formation is apparently unique 
to CeNiSn and CeRhSb; the other Kondo insulators 
SmB & , Ce 3 Bi4Pt 3 and YbB\l display a well-formed gap. 
Curiously, these materials are cubic, leading us to specu- 
late that their higher symmetry prevents the dynamically 
generated contribution to the crystal field from exploring 
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FIG. 2. Normalized thermal conductivity versus tempera- 
ture along the z-axis (solid line) and in the basal plane (dashed 
line). Top: for the Ikeda-Miyake state. Bottom: for the 
quasi-octahedral scenario. Insets show density of states as a 
function of the energy. The adjustable parameters have been 
chosen as V/D = 0.08 and an impurity scattering phase-shift 
of tt/2. 

the region of parameter space where a node can develop. 
At present, we have not included the effect of a mag- 
netic field, which is known to suppress the gap nodes. 6 
There appears to be the interesting possibility that an 
applied field will actually modify the dynamically gener- 
ated crystal field to eliminate the nodes. Finally, we note 
that since the spin-fluctuation spectrum will reflect the 
nodal structure, future neutron scattering experiments 22 
should in principle be able to resolve the axial or octahe- 
dral symmetry of the low energy excitations. 

To conclude, we have proposed a mechanism for the 
dynamical generation of a hybridization gap with nodes 
in the Kondo insulating materials CeNiSn and CeRhSb. 
We have found that Hunds interactions acting on the 
virtual 4/ 2 configurations of the Cerium ions generate 



a Weiss field which acts to co-operatively select a semi- 
metal with nodal anisotropy. Our theory predicts two 
stable states, one axial , the other quasi-octahedral in 
symmetry. The quasi octahedral solution appears to be 
the most promising candidate explanation of the vari- 
ous transport and thermal properties of the narrow-gap 
Kondo insulators. 
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